Magnetic effect from non-magnetic impurity in superconducting CuC>2 plane 
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We propose a new model for impurity center formed by a cathion substitute for Cu, like Zn, in 
Cu02 planes. Its main effect on superconducting electrons is due to the non-zero exchange field 
on O sites, neighbors to the (non-magnetic) impurity. We discuss a strong suppression of d-wave 
order parameter, a zero-energy resonance in local density of states, and spin polarization of charge 
carriers, which can be related to the experimentally observed effects in Zn-doped copper oxides. 
These results are obtained without using the unitary scattering limit. 
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It is recognized that there are two types of impuri- 
ties in high-Tc superconductors (HTSC): i) intrinsic or 
"own", and ii) extrinsic or "foreign". The first type are 
the heterovalent impurities or oxygen vacancies, that is 
the dopants. They supply charge carriers into insulating 
antiferromagnetic (AFM) cuprate planes enabling their 
metallization but also they act as scattering centers 
for carriers. Previously, we showed 0, |l Q{ that forma- 
tion of superconducting (SC) condensate, either of s- and 
d- symmetry, is possible at low enough concentration c of 
such impurities but is prevented by growing fluctuations 
of SC order parameter at higher c. 

The second type are the homovalent impurities, which 
only produce scattering of existing charge carriers and so 
can depress the SC properties of HTSC systems. They 
are well studied in common SC metals, where it was 
stated by Anderson H that non-magnetic impurities have 
practically no effect on SC characteristics. At the same 
time, even low concentration of paramagnetic ions can 
completely destroy SC order, as initially shown in the 
Born approximation by Abrikosov and Gor'kov J(| and 
then confirmed by many authors under more general 
scope § 1 1. 

But in the case of SC copper oxides, an apparent vio- 
lation of this so well theoretically based phenomenolog- 
ical principle was detected. Thus, introducing of non- 
magnetic Zn 2+ ions instead of Cu 2+ into the cuprate 
planes has a suppression effect on HTSC not weaker but 
rather stronger than that by magnetic Ni 2+ ions [TT| . 
This triggered an idea of viewing the non-magnetic impu- 
rity ions in HTSC as extremely strong scatterers | |l2| | so 
that their perturbation potential Vi mp is the biggest en- 
ergy parameter, treated in the unitary limit: Vi mp /W 3> 
1 (W the bandwidth) . This concept was extensively elab- 
orated [l5[ , the principal conclusions being the 
finite density of quasiparticle states (DOS) at the Fermi 
level £p: p(e — > £f) — > p u 7^ 0, and the universal value of 
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quasiparticle conductivity a(ui — > 0) — > a u ^ 0. How- 
ever, apart from the still existing controversies about 
those predictions JlfJ, it should be noted that, unlike 
the dopants, the foreign impurity centers are formed in 
the CuCh plane by a homovalent substitution (as Zn 2+ 
or Ni 2+ for Cu 2+ ), and it is problematic how they could 
produce such a strong perturbation potential. Also we 
notice that the heterovalent non-magnetic scatterers by 
dopants can not produce such effects ||. 

This Letter is aimed on an alternative approach to the 
problem of foreign impurities. It will be shown that ir- 
respectively of the type (magnetic or non-magnetic) of 
the cathion substitute in CuC>2 plane, the resulting cen- 
ter acts on charge carriers as magnetic. In accordance 
with the generally accepted notion, such center should in 
fact strongly suppress SC order either of s- or rf-type, as 
was first qualitatively stated yet by Mahajan et al |fL7|| . 
We note that similar views on the effect of Zn impu- 
rities in HTSC cuprates were expressed in some recent 
publications |fLs| , [l9| , |2p|| , though still focused on unitary 
scattering. Below we consider the problem of isolated 
non-magnetic impurity ion in a CuC>2 plane and its local 
effects on the d-wave SC order parameter, local density of 
states (LDOS), and the itinerant spin polarization. Our 
treatment does not need using the unitary limit, never- 
theless the effects can be quite strong. 

Fig. 1 shows a cathion impurity substitute for Cu in 
a Cu02 plane, like real Zn, Fe, or Ni impurities in high- 
Tc compounds. Associating the charge carriers mostly 
to 0~ holes, we conclude that the main perturbation by 
such impurity (both magnetic and non-magnetic) is due 
to the fact that its neighbor O sites occur in a non-zero 
exchange field by Cu 2+ ions |f2l| , p2| , which is equivalent 
to the effect of magnetic impurity in a common super- 
conductor. On the other hand, there are no reasons to 
consider any sizeable spin-independent perturbation from 
such isovalent impurity. The respective model Hamilto- 
nian is H = H sc + H c + Hi nt . The unperturbed SC term 
H S c = J2k ^k(£kT3 + Ak'TijV'k couples the Nambu spinors 
T>°-k,j.) with Pauli matrices f 4 . The normal 
dispersion £k = W(2 — cos ak x — cos ak y )/A — £f leads to 
the density of states (DOS) po = A/(ttW). The d-wave 
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gap function A k = A6*(ep — £ k )7 k /7m includes the BCS- 
shell factor 6(e^ — £ k ) with the "Debye energy" Ed, the 
symmetry factor 7k = cos ak x — cos ak y with maximum 
absolute value 7,„ = 7reppo, and the gap parameter 

A = ^ 1 E^(a-k,i«k, T ) (1) 

k 

where V is the attraction between two carriers with op- 
posite spins on neighbor O sites. The term H c — ~hS z 
models the (AFM) correlation between the impurity cen- 
ter and its environment, where h ~ Jdd, the Cu-Cu ex- 
change constant, and S is the spin of a fictitious "mag- 
netic impurity". It can be seen as a cluster of four 1/2 
spins of Cu nearest neighbors to real non-magnetic impu- 
rity (Fig. 1). In reality, its quantization axis z is only de- 
fined over time periods no longer than r s ~ h^ s / (a Jdd) ~ 
10 -13 s for spin correlation length £ s ~ a/s/x |2^] and 
doping levels x ~ 0.1 (this also agrees with the NMR 
data ||18|]). However this t s is much longer than typical 
electronic times ~ h/s-p ~ 10 s for HTSC compounds. 
For h > we have (S z ) = s and < s < S, which ac- 
counts for the short-range AFM order, whereas s — ► in 
the paramagnetic limit f3h <C 1. 

We separate the spin-dependent interaction between 
charge carriers and impurity into three parts: 

Hm = H^f + Hf nt + H^ nt , (2) 

where 

k,k' cr=± 

is the "mean-field" (MF) polarization of carrier spins by 
the impurity center, and 

H Lt = JN ^ X! H a jM a jM ,(T ( S * - s)al, a a k:<T , 

k,k' a=± 
k : k' cr=± 

are their interactions with longitudinal and transversal 
fluctuations of S. In the paramagnetic limit: s — > 0, Eq. 
^ is reduced to the common Kondo interaction pi| , p5[ . 
For definiteness, the Cu-0 p-d exchange parameter J is 
considered positive. The functions 

ai,k = 2 cos — cos— , a 2 ,k = 2 cos — sin— , 

x aky , Qjk x , aky 

"3,k = 2 sin — cos— , a 4 , k = 2 sin — sin— , 

realize ID irreducible representations of the C^ v point 
group (of the plaquette surrounding the impurity), so 
that iV -1 J2k a j,k a iLk — fijj'- Distinctive features of the 
perturbation, Eq. |[ compared to the commonly used 
impurity models, are: i) its spatial extension expressed by 



the factors oy,k, ii) additional degrees of freedom by spin 
S, and iii) coupling of S to the local AFM correlations. 

In principle, this impurity center can produce yet an- 
other perturbation, due to a possible role of AFM cor- 
related Cu 2+ spins in the SC coupling between charge 
carriers. Lacking one such spin would locally perturb 
the AkTi term in H sc by some expansions in a^kO^k'. 
This can influence the SC order, alike the simpler case of 
point-like perturbation of s-wave SC coupling [p6[. How- 
ever, for simplicity, we leave this kind of perturbation for 
a separate study. 

We calculate the averages, like Eq. [j], by simple spec- 
tral formula at T = 0: 

(ab) =tt- 1 f F %((b\a)) £ de, (3) 
Jo 

where 37(e) = \im s ^ Q [f(e — iS) — f(e + i5)]/2 and 
((b\ a ))s±iS are the retarded and advanced two-time Green 
functions (GF's). The relevant GF matrix Gk,k' = 
((V'klVv)) m absence of impurity perturbation (J = 0) 
is momentum-diagonal: Gk.k' = ^k.k'Gk, with Gk = 
(e — ^k^3 — Akfi) -1 . The same expression holds for 
the momentum-diagonal GF Gk,k in presence of sin- 
gle impurity, whose effect ~ 1/N is negligible for this 
quantity. However it is only this small impurity effect 
that gives rise to a momentum-non-diagonal GF's Gk,k'- 
They are found from the equation of motion G k ,k' = 
JN' 1 Ek"j *j,k&(s& v , tV +&Q + 6Q,)a jtV , includ- 
ing three scattered GF's: the MF one Gk",k', the longi- 
tudinal G k „ k , = ((ipk"(S z — s)|^! k ,)) and the transversal 

G k -|k' = with 4 = (4 >; ,a_ k)T ). The 

two last terms are analogous to the well known Nagaoka's 
T-term || and treating them with a similar decou- 
pling procedure gives: 

^k*k' = a j,kGkG k ",k'a-, :k ", 
k",j 

^k,k' — X/ a i>kG k (e + /i)A k "Gk»,k'aj,k", 

k"J 

where S 2 = (S 2 z )-s 2 , A k = 5 (5+l)-s(s + l)-E 2 +(l + 
2^k/-^k)T3, Ek = \/£ k +~A k ~, and one ener gy argument 
is shifted: e — ► £ + h, due to the AFM stiffness. Finally, 
we obtain the decoupled equation of motion: 

G k>k ' = N- 1 a jM Gk[Js + J 2 (E 2 G J - + X J -)]G ! k»,k'a j , k » 

k",i 

where G 3 = N' 1 ^ k a 2 k G k , Jtj = N' 1 ^ k a?, k G k (e + 
/i)A"k, and its standard iteration yields in the result: 

G k ,k< = N- 1 J^^kGkfjG^aj, v , (4) 

3 

with the partial T-matrices fj = [Js+J 2 (S 2 Gj+Xj)][1- 
Js — J 2 (S 2 Gj + Xj)] -1 . By the definition of our model, 
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the parameter Js is positive. It is interesting to trace the 
behavior of Tj in the two characteristic limits for AFM 
correlations between Cu 2+ spins. 

In the paramagnetic limit: h — ► 0, s — ► 0, we 
have £ 2 -> 5(5 + l)/3 and X 3 -> 25(5 + l)/3 - 
Ek^.kC 1 + 2Ck/^k)G k f 3 . In neglect of the small 
last term we arrive at: 



tj -► J 2 5(5 + l)Gj[l - J 2 5(5 + l)Gj]- 



generalizing the known results [g, [| for the case of ex- 
tended impurity center. 

Another limit, fully polarized, h oo, s — > 5, corre- 
sponds to S 2 — » 0, Xj — > and results in 

-» J5(l - JSC,)- 1 (5) 

which is only due to the effect of MF magnetic scattering. 
The obvious validity condition for this limit, JS ^> ksT, 
well applies in the SC phase at T < T c ~ A/fc B , so we 
use Eq. [sj for the T-matrices in what follows. 

The local SC correlation is characterized by the aver- 
age A12 = 2V(as ll ias 2 .^) (see Fig. 1) where a site oper- 
ator a n r7 is expressed through band operators: a n , a — 
J^ k e* k ' n ak,(7- Since the phase of Ak is chosen 
zero, A12 is real. For J = 0, this average does not differ 
from the uniform gap parameter: 



A12 



^5>_ k , iak , T )e-(^> 

k 



whereas for J ^ the maximum perturbation of SC order 
near the impurity is given by the suppression parameter 
Vsup = 1 — A12/A. Its value is confined between (for 
pure SC) and 1 (for complete local suppression of SC 
order), and it only results from non-diagonal GF's: 



(«-k,i%M}e !(k - fc - kUl) = (6) 



k.k'^k 



V 
2ttA 



./ — cx 



cfelm TrGj f„- Gj fi , 



where the trace is in Nambu indices and Eqs. ||,H were 
used. Expansion of Gj in Pauli matrices only contains 
fi terms at j — 2,3 (by the parity of <x,,k with respect to 
the permutation k x <-> k y ): JSG2,3 = A + Br? ± Cfi. 
Hence only j = 2, 3 actually contribute in Eq. |6| by 



_ Ve F p 2 f £ - F(e) 

Vsup — A I ' 
4 JO £ 



(7) 



where 

F(e) 



16e 



7tJ5£fA/9, 



2 Im[l 



(1 - .4) 2 - B 2 - C 2 



\C, 



and the complex coefficients A, B, C as functions of en- 
ergy are estimated in the relevant range A < e < £f, 



setting -Bk ~ ~ Wa(k 2 
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JSpoe ef 
JSpo 



C 



B w 

ir J5/?o£fA 
£e 



[In 



fc|)/8: 



:I + 2itt), 

£ F -£ 2 
(I^-ef) 2 ' 



In 



£ D - e)]. 



Numeric analysis of these expressions with realistic pa- 
rameter values: W ~ 2 eV, J5po ~ 1, £f ~ 0.3 eV, 
£d ~ 0.15 eV, shows that the main contribution into the 
integral, Eq. fj], comes from the BCS shell A < e < £d, 
where we have: F(e) « 1, < 1 — F(e) <C 1, while 
< -F(e) < 1 out of this shell, £d < £ < £f (Fig. 
2). Taking in mind the relation ln(e D /A) s» W/ (^7 m ) = 
W 2 /(4y£ F ) which follows from Eq. @, we obtain from 
Eq. |7] almost complete suppression: r\ sup ~ 1. A small 
residual part < 1 — 7? sup <C 1 is only due to a small 
negative deviation of F(e) from unity within the shell 
and to a small negative out-of-shell contribution. Thus, 
for the above indicated choice of parameters we have 
Vsup ~ 96%. The value of 1 — rj sup yet diminishes with 
growing JS, but it should be stressed that no unitary 
limit JSpQ ^ 1 is needed to get such a strong effect. 

The decay of this maximum effect with separation R 
from the impurity is given, in similarity with Eq. g, by 



V 

?7sup(R) = T~ 
7TA 



d£lmTrG2(R)r 2 G 2 (R)f 1 , (8) 



A^Ek^^.kGk, 



and 



Here the matrix Gj(R) 
for R >• £± = ay/W/[8(ep ± e)] it is mainly con- 
tributed by two saddle points in the complex /c-plane: 
±£± —iRr 1 , hence all its matrix elements decay asymp- 
totically like cos(i£/£±) / \J 'R/£±, and the non-diagonal 
elements contain yet the anisotropic factor (A/e)cos2<0 
where ip — arctan R y /R x . However, for the energies 
£ ~ £d relevant here, such anisotropy is less pronounced 
than that in the limit e — > considered by Balatsky 
et al J2?]]. Integration in Eq. || results in asymp- 
totic r) SU p(R) « -q S up \/W/ (8e F ) J2i=± ( u i + fi co ^^ + 
hi cos Aif) (a /R) where hi <C fi ~ Ui ~ 1. This an- 
gular dependence resembles that for LDOS around Zn 
impurity, suggested by Haas and Maki from continuous 
Bogolyubov-de Gennes equations [Q, while the anoma- 
lously slow radial decay should enhance the overall sup- 
pression of SC order. 

Besides the considered local suppression of the 
order parameter, related to the non-diagonal (in 
Nambu indices) elements of GF's Gk.k', there are 
also local effects related to their diagonal elements. 
Thus, the variation of local DOS (LDOS): p(n) = 

(7rA r )" 1 X; kk , #k ImTre l(k_k ' ) ' n Gk,k'T3, attains its max- 
imum at n = 8, nearest neighbor sites to the impu- 
rity, and Eq. || mainly contributes there by j = 1: 
p{5) w ImTrGi(<5)TiGi(<5)f 3 . The relevant GF's are 



4 



GAS) » 2A^ 1 £ k G k = 2(g Q +g 3 T 3 ), where 5o « 
Po^I^f 1 + i(7r/A) arcsin(A/e)] and 33 » poIi^WV^f) 
||. Considered as a function of energy, displays 
a very sharp resonance in the denominator of T\: Re[l — 
2JS , ff0 (e)] 2 ~ (2JSg 3 ) 2 -> at e -» (cf. Fig. 3) 
with the observed peak in the related tunnel conductiv- 
ity j2!|) if the impurity perturbation parameter J is close 
to J cr = l/{2Sg 3 ). This refers to a fine tuned rather 
than unitary perturbation J and agrees with its choice 
made to estimate r) sup . The j = 1 contribution also 
dominates in the Kondo-like local polarization of itin- 
erant spins: m(5) » ^ Mm Tre^^"^ ($)fi<2i($), 
which should explain the observed enhancement of ex- 
change fields on 63 Cu @ and 89 Y JT| nuclei close to Zn 
impurities. A more detailed treatment of these phenom- 
ena will be presented elsewhere. 

The proposed model can be equally applied to isova- 
lent substitutes for Cu, which are magnetic themselves, 
as Ni 2+ or Fe 2+ . But, since the net MF on neighbor 
sites in this case is due to incomplete AFM compen- 
sation of exchange fields by different magnetic ions, the 
perturbation parameter JS may be weaker than that for 
non-magnetic Zn 2+ and so the resulting suppression of 



SC order, as is observed in the experiment 

In conclusion, we developed a microscopic model of 
spin dependent perturbation on charge carriers in CuC-2 
planes, produced by a non-magnetic substitute for Cu. 
An almost complete suppression of d-wave order param- 
eter at nearest neighbor sites to the impurity atom is ob- 
tained, as a result of parallel alignment of carrier spins in 
the exchange field JS by non-compensated Cu 2+ spins, 
and this strong effect is achieved with moderate JS val- 
ues. It decays with distance from impurity rather slowly, 
which can explain the fast destroying of SC order in 
cuprates already at low Zn concentration. The model 
also provides explanation for other local effects, such as a 
sharp resonance of LDOS and local polarization of charge 
carrier spins close to impurity. 
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Figure 1: Effective magnetic perturbation for charge carriers on nearest neighbor 
sites Si to the non-magnetic impurity substitute (grey circle) for Cu 2+ in CuC>2 
plane 




Figure 2: The function F(e) for parameter values W = 2 eV, JSpo = 1, £f = 0.3 
eV, £d = 0.15 eV (the discontinuity is due to the assumed sharp BCS-shell 
boundary). 
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Figure 3: Local DOS on nearest neighbor sites to Zn impurity site, as a 
function of energy p(e), presents a sharp zero-energy resonance at the choice of 
perturbation parameter JSpo = 1.26 ps J cr S. 
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